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The Hund's-rule-exchange induced and coexisting spin-triplet paired and magnetic states are 
considered within the doubly degenerate Hubbard model with interband hybridization. The Hartree- 
Fock approximation combined with the Bardeen-Cooper-SchriefFer (BCS) approach is analyzed for 
the case of square lattice. The calculated phase diagram contains regions of stability of the spin- 
triplet superconducting phase coexisting with either ferromagnetism or antiferromagnetism, as well 
as a pure superconducting phase. The influence of the inter-site hybridization on the stability of the 
considered phases, as well as the temperature dependence of both the magnetic moment and the 
superconducting gaps, are also discussed. Our approach supplements the well known phase diagrams 
containing only magnetic phases with the paired triplet states treated on the same footing. We also 
discuss briefly how to include the spin fluctuations within this model with real space pairing. 
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I. INTRODUCTION 

The spin-triplet superconducting phase is believed to 
appear in Sr2Ru04^, UGe2^, and URhGe^. In the last 
two compounds the considered type of superconducting 
phase occurs as coexisting with ferromagnetism. Addi- 
tionally, even though U atoms in this compounds con- 
tain 5f electrons responsible for magnetism, this multiple- 
band system can be regarded as weakly or moderately 
correlated electron system, particularly at higher pres- 
sure. Originally it has been suggested via a proper quan- 
titative analysis^H^, that the the intra atomic Hund's 
rule exchange can lead in a natural manner to the co- 
existence of superconductivity with magnetic ordering - 
ferromagnetism or antiferromagnetism. 

The coexisting superconducting and magnetic phases 
are discussed in this work within an orbitally degen- 
erate two-band Hubbard model using the Hartree-Fock 
approximation (HF), here combined with the Bardeen- 
Cooper-Schrieffer (BCS) approach, i.e., in the vicinity 
of the ferromagnetism disappearance, where also the su- 
perconductivity occurs. The particular emphasis is put 
on the appearance of superconductivity near the Stoner 
threshold, where the Hartree-Fock-BCS approximation 
can be regarded as realistic. This type of approach can 
be formulated also for other systems^. 

The alternative suggested mechanism for appearance 
of superconductivity in those systems is the pairing me- 
diated by ferromagnetic spin fluctuations, which can 
also appear in the paramagnetic or weakly ferromagnetic 
phase^. Here, the mean- field approximation provides not 
only the starting magnetic phase diagram but also a re- 
lated discussion of the superconducting states treated on 
equal footing. In this approach the spin-fluctuation con- 



tribution appears as a next-order contribution. This is 
the reason for undertaking a revision of the standard 
Hartree-Fock approximation. Namely, we concentrate 
here on the spin triplet states, pure and coexisting with 
either ferromagnetism or antiferromagnetism, depending 
on the relative magnitude of microscopic parameters: the 
Hubbard intra- and inter-orbital interactions, U and U' ^ 
respectively, the Hund's rule ferromagnetic exchange in- 
tegral J, the relative magnitude of hybridization and 
the band filling n. The bare band width W is taken 
as unit of energy. In the concluding Section we discuss 
briefly, how to outline the approach to include also the 
quantum fluctuations around this HF-BCS (saddle point) 
state, as a higher-order contribution. 

The role of exchange interactions is crucial in 
both the so-called t-J model of high temperature 
superconductivit}'^ and in the so-called Kondo-mediated 
pairing in heavy fermion systems^^. In this and the 
following papers we discuss the idea of real space pair- 
ing for the triplet-paired states in the regime of weakly 
correlated particles and include both the inter-band hy- 
bridization and the corresponding Coulomb interactions. 
We think that this relatively simple approach is relevant 
to the mentioned at the beginning ferromagnetic super- 
conductors because of the following reasons. Although 
the effective exchange (Weiss-type) field acts only on the 
spin degrees of freedom, it is important in determining 
the second critical field of ferromagnetic superconduc- 
tor in the so-called Pauli limit as the orbital effects 
in the Cooper-pair breaking process are then negligi- 
ble. The appearance of a stable coexistent ferromagnetic- 
superconducting phase means, that either Pauli limiting 
situation critical field has not been reached in the case 
of spin-singlet pairing or else, the pairing has the spin- 
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triplet nature, without the component with spin S*^ = 0, 
and then the Pauh hmit is not operative. 

The present model with local spin-triplet pairing has 
its precedents of the same type in the case of spin-singlet 
pairing, i.e., the Hubbard model with U < G^^, which 
played the central role in singling out a nontrivial char- 
acter of pairing in real space. Here, the same role is being 
played by the intra-atomic (but inter-orbital) ferromag- 
netic exchange. We believe that this area of research 
unexplored so far in detail opens up new possibilities 
of studies of weakly and moderately correlated magnetic 
superconductor^^. 

The structure of this paper is as follows. In Section 
2 we define the model and the full Hartee-Fock-BCS ap- 
proximation (i. e. mean field approximation for mag- 
netic ordering with the concomitant BCS-type decou- 
pling) for the coexistent two-sublattice antiferromagnetic 
and spin-triplet superconducting phase (cf. also Ap- 
pendix A for details). For completeness, in Appendix 
B we include also the analysis of a simpler coexistant 
superconducting-ferromagnetic phase. In Section 3 we 
provide a detailed numerical analysis and construct the 
full phase diagram on the Hund's rule coupling-band fill- 
ing plane. We describe also the physical properties of the 
coexistent phases. In Appendix C we sketch a system- 
atic approach of going beyond Hartree-Fock approxima- 
tion, i.e., including the spin fiuctuations, starting from 
our Hartree-Fock-BCS state. 



II. MODEL AND COEXISTENT 
ANTIFERROMAGNETIC-SPIN-TRIPLET 
SUPERCONDUCTING PHASE: 
MEAN-FIELD-BCS APPROXIMATION 

We start with the extended orbitally degenerate Hub- 
bard Hamiltonian, which has the form 



il ill'il^l') ^ 



(1) 



where / = 1,2 label the orbitals and the first term de- 
scribes electron hopping between atomic sites i and j. 
For / ^ I' this term represents electron hopping with 
change of the orbital (inter-site, inter-orbital hybridiza- 
tion). Next two terms describe the Coulomb interac- 
tion between electrons on the same atomic site. However 
as one can see the second term contains the contribu- 
tion that originates from the exchange interaction (J). 
The last term expresses the (Hund's rule) ferromagnetic 
exchange between electrons localized on the same site, 
but on different orbitals. This term is regarded as re- 
sponsible for the local spin-triplet pairing in the subse- 
quent discussion. The components of the spin operator 



Sii = {S^i^ S'^i^S^i) used in ^ aquire the form 

-hlia^^y^.hii, (2) 



where (Jx^y^z are the Pauli matrices and hj^ = 
(aj^^,aj^). In our considerations we neglect the 
interaction-induced intra-atomic singlet-pair hopping 
( JaJ-L^aJ-^^ai2|<^i2t + H.c.) and the correlation induced 

hopping {Vnia{a\^a2a + a^^aia) + 1^2 )^, as they 
should not introduce any important new qualitative fea- 
ture in the considered here spin-triplet paired states. 



What is more important, we assume that tlj = tfj and 
t}j = tfj = tij^ i.e., the starting degenerate bands have 
the same width {the extreme degeneracy limit) ^ as we are 
interested in establishing new qualitative features to the 
overall phase diagram, that are introduced by the mag- 
netic pairing. 

As has already been said, the aim of this work is to ex- 
amine the spin-triplet superconductivity coexisting with 
ferromagnetism and antiferromagnetism as well as the 
pure spin-triplet superconducting phase and the pure 
magnetically ordered phases. Labels defining the spin- 
triplet paired phases (A and Al) that are going to be 
used in this work correspond to those defined for su- 
perfiuid ^He according to the Refs.^^ and^. Namely in 
the A phase the superconducting gaps that correspond 
to Cooper pairs with total spin up and down are equal 
(Al = A_i ^ 0, Ao = 0), whereas in the Al phase the 
only nonzero superconducting gap is the one that corre- 
sponds to the Cooper pair with total spin up (Ai ^ 0, 
A_i = Ao = 0). In this section we show the method 
of calculations that is appropriate for the superconduct- 
ing phase coexisting with antiferromagnetism, as well as 
pure superconducting phase of type A and pure antifer- 
romagnetic phase. The corresponding considerations for 
the case of ferromagnetically ordered phases and super- 
conducting phase Al are deferred to the Appendix B. 

From the start we make use of the fact that the full 
exchange term can be represented by the real-space spin- 
triplet pairing operators, in the following manner 



J 



E 

ill'il^l') 



Sir Sij/ 



which are of the form 



it ^ 



TnuUir 



2J^74j^A^^, (3) 



m = 1, 
m = —1, 
0. 



(4) 



Furthermore the inter-orbital Coulomb repulsion term 
can be expressed with the use of spin-triplet pairing op- 
erators and the spin-singlet pairing operators in the fol- 
lowing manner 

([/' + J) ^ = ([/' + J)(^ BjBi + J2 ALAira), 

i i im 

(5) 
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where 



B] 



1 

71 



(6) 



are the inter-orbital, intra-atomic spin-singlet pairing op- 
erators in real space. Using Eq. ([3| and ([5| one can write 
down our model Hamiltonian in the form 



+ [u' + J) ^ b\b. - (J - U') J2 ^L^i 



(7) 



regime, where the metallic state is stable, no orbital or- 
dering is expected (cf. Klejnberg and Spalek in^). 

In our considerations the antiferromagnetic state rep- 
resents the simplest example of the spin- density- wave 
state. In this state, we can divide our system into 
two interpenetrating sublattices A and B. The aver- 
age staggered magnetic moment of electrons on each of 
the N/2 sublattice A sites is equal, < >=< S\ >, 
whereas on the remaining N/2 sublattice B sites we have 



< Sf >=< SI, >E 



< >. In accordance with this 



division into two sublattices, we define different annihi- 
lation operators for each sublattice, namely 



It should be noted here that for J <U' ^ the inter-orbital 
Coulomb repulsion suppresses the spin-triplet pairing 
mechanism and the superconducting phases will not ap- 
pear in the system in the weak-coupling (Hartree-Fock) 
limit. For 3(i electrons^^, U' = U — 2 J, thus the neces- 
sary condition for the pairing to occur in our model is 
U < 3J. Usually, for 3d metals we have U ~ 3 J, so 
it represents a rather stringent condition for the super- 
conductivity to appear in that situation. We use this 
relation to fix the parameters for modeling purposes, not 
limited to 3d systems. This is also because e.g. 5f elec- 
trons in uranium compounds lead to a similar behavior 
as do 3d electrons. One should note that the considered 
here pairing is based on the intra-atomic inter-orbital fer- 
romagnetic Hund's rule exchange. A simple extension to 
the situation with nonlocal J has been considered by X. 
Dai et al.^. Also, as we consider only weakly correlated 



CiiiaA for i G A, 
diiaB for i e B. 



(8) 



The same holds for the creation operators, aj^^. We as- 
sume that the charge ordering is absent. In this situation, 
we can write that 



< riiiA >=< ni2A >=< riiiB >=< ni2B >= n/2, 



(9) 



(10) 

where n is the band filling. In what follows, we treat the 
pairing and the Hubbard parts in the combined mean- 
field-BCS approximation. In effect, we can write down 
the Hamiltonian transformed in reciprocal (k) space in 
the form: 



Hhf -I^N = Y^ 



\{U^2U' -J)-ii 



kll'il^l'a) \^l(T 



hmBJ 



(11) 



k,m=±l 



k,m=±l 



N 



2{J -U') 



where / = + J is the effective magnetic coupling con- 
stant and eki = ek2 = ^k is the dispersion relation. The 
results presented in the next section have been carried 
out for square lattice with nonzero hopping t between 
nearest neighbors only. The corresponding bare disper- 
sion relation in a nonhybridized band acquires the form: 

= —2t cos kx — 2t cos ky . (12) 

As we are considering the doubly degenerate band model 



situation, we make a simplifying assumption that the 
hybridization matrix element ei2k = /^/i^k, where Ph G 
[0, 1] is the parameter, which specifies the hybridization 
strength (i.e. represents a second scale of electron ener- 
gies, in addition to ek). This means that we have just 
one active atom per unit cell with a doubly degenerate 
orbital of the same kind (their spatial asymmetry is dis- 
regarded). One should note that the sums in ([n]) (and 
in all corresponding equations below) is taken over N/2 
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independent k states. In the Hamiltonian written above 
we have also introduced the superconducting spin-triplet 
sublattice gap parameters 



^±1A(B) 



2{J-U') 

N 

2{J-U') 
V2N 



k 

< ^k,oA(B) > • 



(13) 



The terms: 



2U' - J) and f (/7 + 2U' - J) in 
( |ll| ) can be neglected, as they lead only to a shift of the 
reference point of the system energy. One should note 
that since the real-space pairing mechanism is of intra- 
atomic nature, there is no direct conflict with either ferro- 
or antiferro-magnetic ordering coexisting with it. 



A. Antiferromagnetic (Slater) subbands 



The diagonalization of the Hamiltonian (11) can be 
carried out in two steps. In the first step we diagonal- 
ize the one particle part of the Hartree-Fock Hamilto- 
nian (the first two sums of (pT])). Note that we have 
to carry out this step first, since we assume the bands 
are both hybridized and contain pairing part. By in- 
troducing the four-composite fermion operator fj[^ = 

(4iaA' A^2aA^ 42^^)' ^c cau cxprcss the ouc par- 

ticle Hamiltonian in the following form 



(14) 



where fk 



{^1)1 and 
/ 



V 



-cjISl 

ek 

ekl2 



ka 



_ 

-cjISl 
eki2 
ek 



ek 
eki2 





eki2 \ 

ek 

0_ 
cjISl j 



(15) 



To diagonalize this Hamiltonian we introduce a general- 
ized Bogoliubov transformation to new operators a]^i(jA 
and a}ziaB in the following manner 



/ <^klcrA \ 
<^k2crA 

<^klcrB 
\ <^k2crB 



-V, 



kcr 



kcr 



k(T 





/ ^klcrA \ 




^k2crA 




^klcrB 1 


/ 


\ ^k2crB / 



where 



(±) 



kcr 



V, 



(±) 



kcr 



1 

1 
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V(ek ± eki2)2 + )2 



1/2 



1/2 



(16) 



(17) 



One should note that the symbols A and B that appear 
as indexes of the new quasi-particle operators di^iaA and 



o^MaB-, single out the new, hybridized, quasi-particle sub- 
bands and do not correspond to the sublattices indices A 
and B as in the case of operators ak/crA and ai^icrB- The 
dispersion relations in the new quasi-particle representa- 
tion acquire the form 



CklA 
CklB 

ek2A 

ek2S 



(ek + eki2)^ + (/5|)2. 



(ek + eki2)^ + (/5|)^ 



(ek-ekl2)^ + (/5|)^ 



(ek-eki2)^ + (/5|)^. 



(18) 



As one can see, the new dispersion relations do not de- 
pend on the spin quantum numbers of the quasi-particle. 
In general if eki2 is not ~ ek, we may have four non- 
degenerate Slater subbands, which is not the case con- 
sidered here. To express the pairing operators that are 
present in the Hamiltonian ( 11 ) in terms of the new quasi- 
particle operators, one can use relations (16) and the 
definitions (|4|. The explicit form of the original pair- 
ing operators in terms of the newly created quasi-particle 
operators is provided in Appendix A. 



B. Quasiparticle states for the coexistent 
antiferromagnetic and superconducting phase 



In the second step of the diagonalization of (11), a 
generalized Nambu-Bogolubov-De Gennes scheme is in- 
troduced to write down the complete H-F Hamiltonian 
again in the matrix form, which allows for an easy deter- 
mination of its eigenvalues. With the help of composite 



creation operator fj[^ = (a) 



4ic 



B^' 



-k2crB), 



we can construct this new 4x4 Hamiltonian matrix and 
write 



HHF-f^N = ^fi[,Hkafka+2^(gk2A+ek2s)-2/i7V+C, 

(19) 



kcr 



with 



Hkcr = 



/ CklA - 

^Ikcr 



61 



^3kcr 



^^Ikcr 

^4kcr 





M ^Ika U ()3ka \ 

1 

J 

(20) 

and fk = (fk)^- The parameters (5/ka are defined as fol- 
lows 





(^4kcr 
CklB — 
^2kcr 



^^3kcr 


^^2kcr 

-ek2B + M 



(5ik. = A^aULU^^ + A^B^k^'.^k.. 



(21) 



^2k. = KaV^^V^, + t/+ , 

(^3k. = -A^aU^V-^ + A^bU-^V^^, 
(54k. = -A^aV^^U^^ + A^bV^-U^^. 

Constant C contains the last two terms of the r. h. s. of 



expression (11). Hamiltonian (19) and matrix (20) have 
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been written under the assumption that Aqa = Aqb = 0. 
Calculations for the more general case of nonzero gap pa- 
rameters for m = have been also done, but no stable 
coexisting superconducting and antiferromagnetic solu- 
tions have been found numerically. The only coexist- 
ing solutions that have been found, fulfill the relation 
Aqa = Aqb = 0. This fact can be understood by the fol- 
lowing argument. As in the antiferromagnetic state all 
lattice sites have nonzero magnetic moment, the Cooper 
pairs in the spin-triplet state for m = (i.e. with the 
total spin = 0, corresponding < A\^o >) are not 
likely to appear. Nevertheless, we present the matrix 
form of the Hamiltonian ( pTj ) for the mentioned most 
general case, in Appendix A. In our considerations here, 
we limit also to the situation with the real gap param- 
eters A^-j^^^^^j = A^i^^^). Then, the straightforward 



diagonalization of (20) yields to the following Hamilto- 
nian 

kZcr 

k 



kcr 



(22) 



where AkZcrA(B) are the eigenvalues of the matrix (20) 
and a\^iaA{B) (<^LcrA(B)) quasi-particle annihila- 

tion (creation) operators, related to the original annihi- 
lation and creation operators a\^icr^ ^k^^r fro^ the first 
step of our diagonalization, via generalized Bogoliubov 
transformation of the form 



fkcr = Uj^^gkcr, 



(23) 



with g]^^ = (aki^^,a-k2aA,Q^kio^<^-k2aB). Eigenvec- 
tors of the Hamiltonian matrix (20) are the columns of 



the diagonalization matrix Uj^. Using the definitions of 
gap parameters A±i^, A±i^, the average number of 
particles per atomic site n = < fiu^A + t^u^a >, 
and the average magnetic moment per band per site 
=< fiii^A — ^iiiA > /2, we can construct the set of 
self-consistent equations for the mean-field parameters 
{A±iA, A±iB, S^) and for the chemical potential. The 
averages that appear in the set of self-consistent equa- 
tions, < <^k^crA(B)^k^crA(B) >7 can be replaced by the cor- 
responding Fermi distribution functions 



/((-1)'+'A: 



k^crA(B) 



) = l/[exp(/3(-l)'+^A: 



^k^crA(B)) + 1], 

(24) 

where /3 = l/ksT. The eigenvalues and the eigenvec- 
tors of ( 20 ) are evaluated numerically while executing the 



numerical procedure of solving the set of self-consistent 
equations. For a given set of microscopic parameters n, 
J, [/, and temperature T, the set of self-consistent 
equations has several solutions that correspond to differ- 
ent phases. Free energy can be evaluated for each of the 



solutions that have been found and the one that corre- 
sponds to the lowest value of the free energy is regarded 
as the stable phase. The expression for the free energy 
functional in the considered case has the form 



F = 



-y 



ln( l+exp(-/?(-l)'+iAk,aA) 



kJcr ^ 

ln( l+exp(-/?(-l)'+iAk;aB) 



(25) 



+ 2 ^(ek2A + ^\^2b) + y^(Ak2crA 

k kcr 

+ Ak2aB)-/i(2-n)A/' + C. 

Numerical results are carried out for square lattice with 
nonzero hopping t between the nearest neighbors only. 
The described above numerical scheme is executed for 
the following selection of phases: 

• normal state (NS): A±ia{b) = 0, = 

• pure superconducting phase type A (A): 
A±iA(B) = A 7^ 0, 5j = 

• pure antiferromagnetic phase (AF): A±i^(^) = 0, 

• coexistent superconducting and antiferromagnetic 
phase (SC+AF): A^^a{b) + 0, SI ^ 

The ferromagnetically ordered phases, that will also be 
included in our considerations in the following Sections, 
are listed below: 

• pure saturated ferromagnetic phase (SFM): 



^±1A{B) 



= 0, 5^ = S:^ 



u{max) 



7^0 



• pure nonsaturated ferromagnetic phase (FM): 
^±iMB) = 0, < 5^ < 

• saturated ferromagnetic phase coexistent with su- 
perconductivity of type Al (Al+SFM): 



AiA(B) = Al 7^ 0, A_i^(B) = 0, 55 = S[ 



{max) 



7^0 



• nonsaturated ferromagnetic phase coexistent with 
superconductivity of type Al (Al+FM): 
^iA(B) ^ Al ^ 0, A_i^(B) = 0, < 5S < 

It should be noted that S^^ refers to the uniform mag- 
netic moment per band, per site in the ferromagnetically 
ordered phases, whereas is the staggered magnetic mo- 
ment that corresponds to the antiferromagnetic phases. 
One could also consider the so called superconducting 
phase of type B for which all superconducting gaps (in- 
cluding AoA(B)) are equal and different from zero. How- 
ever this phase never coexists with magnetic ordering. 
What is more important in the absence of magnetic or- 
dering the superconducting phase A has always lower free 
energy than the B phase. Therefore the superconducting 
B phase is absent in the following discussion. 
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III. RESULTS AND DISCUSSION 



calculated gap parameters fulfill the relations 



We assume that U' = U — 2J and U = 2.2 J, so there 
are actually two independent parameters in the consid- 
ered model - n and J. The energies have been normalized 
to the bare band-width W = 8|t|, and T expresses the 
reduced temperature, T = ksT/W. 



A. Overall phase diagram: coexistent 
magnetic- paired states 

In Fig. [l] a-d we present the complete phase diagrams 
in coordinates (n, J) for different values of the hybridiza- 
tion parameter /3h- They comprise sizeable regions of 
stable spin-triplet superconducting phase coexisting with 
either ferromagnetism or antiferromagnetism, as well as 
pure superconducting phase A. In the phase SC+AF the 



(a) 



(c) 



AF 

A1+SFM 1 A1+SFM 




1 2 3 

BAND FILLING, n 

AF 



0.25 



I 0.2 
_i 

§0.15 
I 0.1 

X 

0.05 



A1+SFM V 


^ A1+SFM 






A1+FM^ 


1(^ai+fJ 






- aAaf 


AF / X I 


'_ NS \ 


/, :n 


^ . 1 1 \ 





A+iA = A_iB = A+, 
A_i^ = A+iB = A_, 
A+ > A_. 



(26) 



For the singlet paired state one would have A+i^ = 
— A_iA, which is not the case here. For the case of 
half filled band, n = 2, the superconducting gaps 
and A_ vanish and only pure (Slater type) AF survives. 
The appearance of the AF state for n = 2 corresponds 
to the fact that the bare Fermi-surface topology has a 
rectangular structure with Q = (tt, tt) nesting. This fea- 
ture survives also for (3^ ^ 0. Also, the symmetry of the 
phase diagrams with respect to half-filled band situation 
is a manifestation of the particle-hole symmetry, since 
the bare density of states is symmetric with respect to 
the middle point of the band. 
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FIG. 1. (Color online) Phase diagrams in space (n, J) for T — and for different values of the fih parameter: (a) fih — 0.00, 
(b) = 0.04, (c) = 0.06, (d) — 0.11. Labels representing different phases are described in main text. One sees that 
practically all magnetic phases here are in fact the coexistent phases with superconductivity except the half filled situation 
where we have pure AF phase 



This feature of the problem provides an additional test 
for the correctness of the numerical results. It is clearly 
seen from the presented figures that the influence of hy- 
bridization is significant quantitatively when it comes to 



the superconducting phase A, as the region of its stability 
narrows down rapidly with the increase of {3^. 

The stability areas of Al+FM and NS phases expand 
on the expense of A and Al+SFM phases. With the 
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further increase of the hybridization, the stabihty of A 
phase is completely suppressed, as shown in Fig. [l]d. The 
regions of stable antiferromagnetically ordered phase do 
not alter significantly with the increasing hybridization. 
To relate the appearance of superconductivity with the 
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FIG. 2. (Color online) Phase diagrams zoomed in space (n, J) 
with the dashed line marking the Stoner threshold for the 
onset of pure ferromagnetism. The values of the hybridization 
parameter are Ph' (a) /Sh — 0.00, (b) I3h — 0.11, while the 
temperature is T = 10~^. 

onset of ferromagnetism we have marked explicitly in Fig. 
[2] the Stoner threshold on the phase diagram. One sees 
clearly that only the Al phase appearance is related to 
the onset of ferromagnetism. What is more important, 
the FM phase coexisting with the paired Al phase, be- 
comes stable for slightly lower J values than the Stoner 
threshold for appearance of pure FM phase. The Al+FM 
coexistence near the Stoner threshold can be analyzed 
by showing explicitly the magnetization and supercon- 
ducting gap evolution with increasing J. This is shown 
m Figs. [3]and|4| One sees explicitly that the nonzero 
magnetization appears slightly below the Stoner thresh- 
old and is thus induced by the onset of Al paired state. 
In other words, superconductivity enhances magnetism. 
But opposite is also true, i.e., the gap increases rapidly 
in this regime, where magnetization changes. The situ- 
ation is preserved for nonzero hybridization. The tran- 
sition A^ Al+FM is sharp, as detailed free-energy plot 
shows. For j3h = 0.11 in a certain range of J the su- 
perconducting solutions Al+FM and A cannot be found 



by the numerical procedure. That is why the curves 
representing the gap parameters A and free energy sud- 
denly break. The most important and surprising conclu- 
sion is that in Al+FM phase only the electrons in spin- 
majority subband are paired. This conclusion may have 
important practical consequences for spin filtering across 
NS/ Al+FM interface, as discussed at the end. Neverthe- 
less, one should note that the partially polarized (FM) 
state appears only in a narrow window of J values near 
the Stoner threshold, at least for the selected density of 
states. 
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FIG. 3. (Color online) Magnetic moment (per orbital per 
site) , ground state energy and superconducting gap as a func- 
tion of J near the Stoner threshold for n — 1 and = 0.0. 
Black vertical line in the inset marks the onset of saturated 
magnetism at the Stoner threshold. 



Summarizing, we have supplemented the well known 
magnetic phase diagrams with the appropriate stable and 
spin-triplet paired states. A relatively weak hybridiza- 
tion of band states destabilizes pure paired states but sta- 
bilizes coexistent superconducting-magnetic phases ex- 
cept for the half- filled band case, when the appearance of 
the Slater gap at the Fermi level excludes any supercon- 
ducting state. A very interesting phenomenon of pairing 
for one-spin (majority) electrons occurs near the Stoner 
threshold for the onset of FM phase and extends to the 
regime slightly below threshold. 
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FIG. 4. (Color online) Magnetic moment (per orbital per 
site) , ground state energy and superconducting gap as a func- 
tion of J near the S toner threshold for n = 1 and /3h = 0.11. 
Black vertical line in the inset marks the onset of magnetism 
at the St oner threshold. 



narrow range of n, a pure AF phase is stable. The inset 
in Fig. [6^ shows that there is aweak first-order transi- 
tion between the AF+SC and SC phases as a function of 
doping. The A phase is not stable in this case. 
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B. Detailed physical properties 

In Figs. [5] and [6] we show the low-temperature values 
of superconducting gaps and the staggered magnetic mo- 
ment as a function of band filling. In the SC+AF phase 
both gap parameters A+ and A_ decrease continuously 
to zero as the system approaches the half filling. On 
the contrary, the staggered magnetic moment Sg reaches 
then the maximum. For the case of (3^ = 0.0, below the 
critical value of band filling, ric ^ 1.45 , the gap param- 
eters and A_ are equal and the staggered magnetic 
moment vanishes. In this regime the superconducting 
phase of type A is the stable one. For A phase the super- 
conducting gap decreases with the band-filling decrease 
and becomes zero for some particular value of n. Below 
that value the NS (paramagnetic state) is stable. It is 
clearly seen that the appearance of two gap parameters 
above ric is connected with the onset of the staggered- 
moment structure, as above Uc we have ^ (cf. Fig. 
^p). For comparison, we also show the staggered moment 
for pure AF in Figs. [sJd and [gJ) (dashed line). As one 
can see, the appearance of SC increases slightly the stag- 
gered moment in SC+AF phase. For = 0.11 below 
some critical value of band filling Uc ~ 1.473 in a very 



FIG. 5. (Color online) Low temperature values of the super- 
conducting gaps and the staggered magnetic moment both as 
a function of band filling for = and J = 0.175. The sta- 
ble phases are appropriately labelled in the regimes of their 
stability. Note that A_ << A+, i.e., the paired state is closer 
to Al state than to A state in the coexistent regime. 

One should mention that the easiness, with which the 
superconducting triplet state is accommodated within 
the antiferromagnetic phase stems from the fact that the 
SC gaps have an intra-atomic origin and the correspond- 
ing spins have then the tendency to be parallel. There- 
fore, the pairs respect the Hund's rule and do not disturb 
largely the staggered- moment structure, which is of in- 
teratomic character. 

In Fig. [7|we show temperature dependence of the free 
energy for the six considered phases, for the set of mi- 
croscopic parameters selected to make the SC+AF phase 
stable at T = and for (3^ = 0. Because the free-energy 
values of the A and NS phases are very close, we exhibit 
their temperature dependences zoomed in Fig. [7]3. The 
same is done for the free energy of phases Al+FM and 
FM. For the same values of n, J, U and [/^ the temper- 
ature dependence of the superconducting gaps and the 
staggered magnetic moment in SC+AF phase are shown 
in Fig. [Sj for selected f3h values. For given f3h below 
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FIG. 6. (Color online) Low temperature values of the su- 
perconducting gaps and the staggered magnetic moment as 
a function of band filling for 13 h = 0.11 and J = 0.175. 
Note the disappearance of the pure A phase and that again 
A- << A+. The inset in (b) illustrates the fact that, a pure 
AF phase appears in a very narrow regime of n before the 
SC+AF phase becomes stable, whereas the inset in (a) shows 
the free energy of those two phases for n close to nc when a 
weak first-order transition occurs. 



Slater gap at Tn makes the superconducting transition 
to disappear. As one can see form Figs. [S] and [9) with 
the increase of the critical temperature Ts is decreas- 
ing slightly while the Neel temperature increases, but the 
ratio remains almost fixed, Tn/Tq ^ 10. 
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FIG. 7. Color online) (a) Temperature dependence of the 
free energy for considered phases, for n = 1.9 and J = 0.175 
when the SC+AF phase is stable at T = 0. The free-energy 
values of A and NS phases are very close, so we exhibit their 
temperature dependence blown up in part (b). 



the superconducting critical temperature T^, the stag- 
gered magnetic moment and the superconducting gaps 
have all nonzero values which means, that we are dealing 
with the coexistence of superconductivity and antiferro- 
magnetism in this range of temperatures. Both and 
A_ vanish at T^, while the staggered magnetic moment 
vanishes at the Neel temperature, T/v » Ts. In Fig. 
[9] one can observe that there are two typical mean- field 
discontinuities in the specific-heat at Ts and Tn for a 
given Ph. The first of them, at T^, corresponds to the 
phase transition from the SC+AF phase to the pure AF 
phase, while the second, at T^v, corresponds to the tran- 
sition from the AF phase to the NS phase. The values 
of the ratios of the specific heat jump (Ac/cn) at Tc 
that correspond to Ph = 0.0, 0.4, 0.6, are 15.075, 16.298, 
17.375 respectively. No antiferromagnetic gap is created, 
since we have number of electrons n < 2. The specific 
heat discontinuity at AF transition is due to the change 
of spin entropy near T^. For n = 2 the formation of the 



Temperature dependence of free energies of relevant 

0) for the micro- 



10 iPh 



phases are presented in Fig. 
scopic parameters selected to make the Al+FM phase 
stable at T = 0. Free energies for A and Al+FM phases 
are drawn only in the low-T regime (Fig. [ToJd) for the 
sake of clarity. The corresponding temperature depen- 
dence of the superconducting gaps, magnetic moment, 
and specific heat in Al+FM phase for three selected val- 
ues of ph are shown in Fig. 11 Analogously as in the 



SC+AF case, the system undergoes two phase transi- 
tions. The influence of hybridization on the tempera- 
ture dependences is also similar to that in the case of co- 
existence of superconductivity with antiferromagnetism. 
With the increasing the critical temperature- is 
decreasing slightly, whereas the Curie temperature, Tc 
is slightly increasing, but still Tc /Ts ~ 5. The values of 
the ratios of the specific heat jump (Ac/cn) at Tc that 
correspond to f3h = 0.0, 0.2, 0.4, are 1.329, 1.421, 0.793 
respectively. 
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tion decreases T^. It should be noted that the values of 
are very close to zero. This is necessary to assume 
for the A phase to appear. The values of the ratios of 
the specific heat jump (Ac/cn) at Tc that correspond 
to f3h = 0.0, 0.035, 0.006, are 1.382, 1.326, 1.202 respec- 
tively. 

In Table [T| we have assembled the exemplary values of 
mean field parameters, chemical potential, as well as free 
energy for two different sets of values of microscopic pa- 
rameters corresponding to the low-temperature stability 
of two considered here superconducting phases: SC+AF 
and Al-I-FM. For the two sets of values of n and J, the 
free energy difference between the stable and first unsta- 
ble phases is of order 10~^. The values for the stable 
phases are underlined. 
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FIG. 8. (Color online) Temperature dependences of the su- 
perconducting gaps A+, A- and of the staggered magnetic 
moment for n = 1.9, J = 0.175 and for selected values of the 
I3h parameter. Note that Ts << Tn- 
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FIG. 9. (Color online) Temperature dependence of the specific 
heat for n = 1.9, J = 0.175 and for selected values of /3h 
parameter. The behavior is almost independent of /3h value 
and the ratio Tn /Ts ~ 10. 
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FIG. 10. (Color online) Temperature dependence of the free 
energy for n = 1.0 and J — 0.31625 when the Al+FM phase 
is stable at T = 0. AF phases does not appear in this case. 
Free energies for A and Al+FM phases are shown in the low- 
T regime (b) for the sake of clarity. 



For the sake of completeness, in Fig. [12] we provide 
the temperature dependence of superconducting gap for 
the values of parameters that correspond to stable pure 
superconducting phase of type A at T = and for three 
different values of ph- In this case, neither the antiferro- 
magnetically ordered nor the pure ferromagnetic phases 
do exist. As in previous cases, the increasing hybridiza- 



IV. CONCLUSIONS AND OUTLOOK 

We have carried out the Hartree-Fock-BCS analysis 
of the hybridized two-band Hubbard model with the 
Hund's-rule induced magnetism and spin-triplet pair- 
ing. We have determined the regions of stability of the 
spin-triplet paired phases with Aq = 0, coexisting with 
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FIG. 11. (Color online) Temperature dependence of the su- 
perconducting gaps A+, A_, magnetic moment (a), and spe- 
cific heat (b), for n = 1.0, J = 0.31625 and for selected values 
of Ph- Qualitative features do not alter appreciably even for 
/3h = 0.4. The ratio Tc /Ts ^ 5. 



TABLE I. Exemplary values of the mean field parameters, 
the chemical potential, and the free energy of the considered 
phases at T = 10~^, for two different sets of values of micro- 
scopic parameters: n, J. The underlined values correspond 
the stable phases. The numerical accuracy is better than the 
last digit. 
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FIG. 12. (Color online) Temperature dependences of the su- 
perconducting gaps A+, A- (a), and the specific heat (b) for 
n — 1.25, J — 0.175 and for selected values of I3h parameter. 



either ferromagnetism (Al+FM) or antiferromagnetism 
(SC+AF), as well as pure paired phase (A). We have an- 
alyzed in detail the effect of inter band hybridization on 
stability of the those phases. The hybridization reduces 
significantly the stability regime of the superconducting 
phase A, mainly in favor of the paramagnetic (normal) 
phase, NS. For large enough value of {Ph > 0.08), 



the A phase disappears altogether. When it comes to 
magnetism, with the increase of the stability regime 
of the saturated ferromagnetically ordered phase is re- 
duced in favor of the non-saturated. The influence of 
the hybridization on the low-temperature stability of the 
SC+AF phase is not significant. When the system is 
close to the half filling, the SC+AF phase is the sta- 
ble one. However, for the half filled band case (n = 2), 
the superconductivity disappears and only pure antifer- 
romagnetic state survives, since the nesting effect of the 
two-dimentional band structure prevails then. 

We have also examined the temperature dependence 
of the order parameters and the specific heat. For 
both coexistent superconducting and magnetically or- 
dered phases (SC+AF and Al+FM) one observes two 
separate phase transitions with the increasing tempera- 
ture. The first of them, at substantially lower temper- 
ature (T^), is the transition from the superconducting- 
magnetic coexistent phase to the pure magnetic phase 
and the second, occuring at much higher temperature 
(T/v or Tc), is from the magnetic to the paramagnetic 
phase (NS). The hybridization has a negative influence 
on the spin-triplet superconductivity, since it reduces the 
critical temperature for each type of the spin-triplet su- 
perconducting phase considered here. On the other hand. 
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the Curie (Tc) and the Neel (Tn) temperatures are in- 
creasing with the increase of the parameter, as it gen- 
erahy increases the density of states at the Fermi level 
(for appropriate band fillings). 

One sholud note that since the pairing is intra-atomic 
in nature the spin-triplet gaps are of the s type. This 
constitutes one of the differences with the corresponding 
situation for superfluid ^He^ where they are of p type^. 

It is also important to note that the paired state ap- 
pears both below and above the Stoner threshold for the 
onset of ferromagnetism (cf. Fig. |2|, though its nature 
changes (A and Al states, respectively). In the ferro- 
magnetically ordered phase only the spin- majority carri- 
ers are paired. This is not the case for AF+SC phase. It 
would be very interesting to try to detect such highly un- 
conventional SC phase. In particular the Andreev reflec- 
tion and in general, the NS/SC conductance spectroscopy 
will have an unusual character. We should see progress 
along this line of research soon. 

As mentioned before all the results presented in the 
previous section has been obtained assuming that U = 
2. 2 J and U' = U — 2J . Having said that the value of 
= J — U' determines the strength of the pairing 
mechanism while / = + J is the effective magnetic 
coupling constant, one can roughly predict how will the 
change in the relations between [/, U' ^ and J result. It 
seems reasonable to say that the larger is J with respect 
to U' the stronger the superconducting gap in the paired 
phases. This would also result in the increase of Tc and 
a corresponding enlargement of the area occupied by the 
superconducting phases on the diagrams. Furthermore 
the increase of U with respect to J should result in the 
increase of the ratios Tc/Ts and Tn/Ts- This is because 
in that manner we make the magnetic coupling stronger 
with respect to the pairing. If we however increase U but 
do not change J^, then the strength of the pairing would 
be the same but the magnetic coupling constant would be 
stronger so this would favor the coexistent magnetic and 
superconducting phases with respect to the pure super- 
conducting phase. Quite stringent necessary condition 
for the pairing to appear J > U' (equivalent to 3J > /7 
if we assume U' = — 2 J, as has been done here) in- 
dicates that only in specific materials one would expect 
for the Hund's rule to create the superconducting phase. 
This may explain why only in very few compounds the 



coexistent ferromagnetic and superconducting phase has 
been indeed observed. Obviously, one still has to add the 
paramagnon pairing (cf. Appendix C). 

It should be noted that more exotic magnetic phases 
may appear in the two-band model^^. Here, we neglect 
those phases because of two reasons. First the lattice 
selected for analysis is bipartite, with strong nesting (AF 
tendency). Second the additional ferrimagnetic, spiral, 
etc., phase might appear if we assumed that the second 
hopping integral t' 0. Inclusion of t' would require a 
separate analysis, as the lattice becomes frustrated then. 
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APPENDIX A. HAMILTONIAN MATRIX FORM 
IN THE COEXISTENT SC-hAF PHASE AND 
QUASIPARTICLE OPERATORS 

In this Appendix we show the general form of the 
Hamiltonian matrix Hk and the pairing operators ex- 
pressed in terms of the quasi-particle creation operators 
from the first step of the diagonalization procedure diss- 
cused in section 2. 

For the case of nonzero gap parameters Aoa(b) we have 
to use eight element composite creation operator 

(^kltA'^kliA'^-k2tA,^-k2iA,ttki^^, 

to write down the Hamiltonian ( pTj ) in the matrix form 

Hhf -^N = Y^ f^Hkfk + 2 ^ (ek2A + e^2B ) - 2/i7V + C, 
k k 

(27) 

where fk = (fk)^ i ^ind 
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The Sikcra-' are the generahzation of parameters intro- 
duced earher in Eq. (16). 



(29) 



where A 



Below we present the pairing operators expressed in 
terms of the quasi-particle creation operators that we 
have introduced during the first step of the diagonahza- 
tion procedure of the Hamiltonian ( [TT] ). 



^kcrA — ^ka^kcr^klcrA^-k2crA + ^k^K:cr^klcrB^-k2aB 
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tonian into the reciprocal space to get 
Hhf - fJ^N = X](ek - M - crIS^)nMa 
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(32) 



where is the uniform average magnetic moment and 
this time the sums are taken over all N independent k 
points, as here we do not need to perform the division 
into two sublattices. In the equation above we have omit- 
ted the terms that only lead to the shift of the reference 
energy. Next, we diagonalize the one particle part of the 
H-F Hamiltonian by introducing quaziparticle operators 



<^klo 



1 



V2 



(ttklo 



<^k2o 



V2 



( — «klo 



<^k2cr), 
<^k2cr), 



(33) 



with dispersion relations 



^kOA — ^/^ / .(^kcr^ka^klcrA^-k2aA 



'^k^'^ko-^klcrB^-k2aB "^kcr ^ka^klcrB^-k2aA 



^ko-K:a^klcrA^-k2aB)' 
^kOB ~ y^C^k^'^ka^klcrA^-k2aA 



^ko-^ka^klcrB^-k2aB + ^kcr^ka^^klcrB^-k2aA 
^k^^ka^klcrA^-k2aB)- 



(31) 



APPENDIX B. HAMILTONIAN MATRIX AND 
QUASIPARTICLE STATES FOR THE 

COEXISTENT 
FERROMAGNETIC-SPIN-TRIPLET 
SUPERCONDUCTING PHASE 



In this Appendix we show briefly the approach to 
the coexistent ferromagnetic-spin-triplet superconduct- 
ing phase within the mean-fleld-BCS approximation. 
Inanalogy to the situation considered in Section 2, we 
make use of relations ([3|, ([5| and transform our Hamil- 



Cklcr 

ek2cr 



ek 



alS' 



|ei2k| 
kl2k| 



(34) 



Using the 4-component composite creation operator 

^k = (^kif ^kii' ^-k2t5 ^-k2i), wc cau coustruct the 4x4 
Hamiltonian matrix and write it in the following form 



Hhf -tJ.N = J2 fk^kfk + Yl ^1^2. + C, 



(35) 



kcr 



where 



Hk = 



/ 


^klt 





Ai 


Ao \ 







Ckii 


Ao 


A_i 






A* 


— ek2t 





v 


A* 
^0 


All 





-^k2i / 



(36) 



with fk = (fk) • Symbol C refers to the last two terms 
of r. h. s. of expression (32). After making the diago- 



nalization transformation of (36) we can write the H-F 
Hamiltonian as follows 



Hhf - IJiN 



\<ilo- 



kcr 



Ak2cr) + 

(37) 

where we have again introduced the quasiparticle opera- 
tors ak^cr and o^k^^. Assuming that Aq = and that the 
remaining gap parameters are real, we can write down 
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the dispersion relations for the quasi-particles AkZa in the 
fohowing way 





- IS^ f 


+ A2 


+ /3/i|ek 


(ek - + 






+ ;S^|ek 




- IS^ Y 


+ A2 




(ek - + 


is^Y + 




-/^^kk 



(38) 



In this manner we have obtained the fuhy diagonahzed 
Hamiltonian analyticahy for the case of superconductiv- 
ity coexisting with ferromagnetism. Next, in the similar 
way as for the antiferromagnetically ordered phases, we 
can construct the set of self consistent equations for the 
mean field parameters A±i, and for the chemical po- 
tential, as well as construct the expression for the free 
energy. 



APPENDIX C. BEYOND THE HARTREE-FOCK 
APPROXIMATION: 
HUBBARD-STRATONOVICH 
TRANSFORMATION 

In outlining the systematic approach going beyond the 
Hartree-Fock approximation we start with Hamiltonian 
^ with the singlet pairing part {U' ^ J)Y,. B\Bi 
neglected, i.e. 

H = Ho^uY^ hutfiiii -J^'Yl AL^im, (39) 



where contains the hopping term, and = J — 
U' . We use the spin-rotationally invariant form of the 
Hubbard term 



(40) 



where hu = hucr and j~Li is an arbitrary unit vector es- 
tablishing local spin quantization axis. One should note 
that, strictly speaking we have to make the Hubbard- 
Stratonovich transformation twice, for each of the last 



two terms in (40) separately. The last term will be effec- 



tively transformed in the following manner 

im 

- ^(iLA,™ + i^A*^ - 

im 

(41) 

where A^^ is the classical (Bose) field in the coherent- 
state representation. The term (|40]) can be represented 
in the standard form through the Poisson integral 



exp - 



1 



dxi( - — 
27r./_^ V 2 



(42) 



In effect, the partition function for the Hamiltonian (39) 
will have the form in the coherent-state representation 

X exp j - / dr\ 4a 

-E 



''ij 



d_ 

a7 



ijW a 

Ai™(r)iL(T) + A*^(T)ii™(r) 



|A,^(r)p 
J" 



(43) 



where we have included only the spin and the pairing 
fiuct nations. In the present paper t-^- = Ujdii' + (1 — 
Sii')tjj. Also the integration takes place in imaginary 
time domain and the creation and anihilation operators 
are now Grassman variables^^. In this formulation A^^ 
and Xi represent local fields which can be regarded as 
mean (Hartree-Fock) fields with Gaussian fiuctuations. 

With the help of (|43]) we can define "time dependent" 
effective Hamiltonian. 



jy(r)^ ^t«:aL(r)a,,v(r)-J^^ 

ijlUa im 



(44) 



where now the fiuctuating dimensionless fields are defined 
as 



Xii[r) = — , /\im[T) 



U 



(45) 



Note that the magnetic molecular field ~ UXii{r) is sub- 
stantially stronger than the pairing field ~ J^Ai^(r). In 
the saddle point approximation Xii{r) = A^^e^, Ai^(r) = 
A*^(r) = A, and we obtain the Hartree-Fock-type ap- 
proximation. Therefore, the quantum fiuctuations are 
described by the terms 

-UY,^Xi{T)-Su{r) 



(46) 



The first term represents the quantum spin fiuctuations 
of the amplitude SXii{r) = Xii{r) — Ae^, the second de- 
cribes pairing fiuctuations. Both fiuctuations are Gaus- 

^ 2 

sian due to the presence of the terms ~ JA^^(r) and 
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\6Aim{r)\'^ . In other words, they represent the higher- 
order contributions and wih be treated in detail elswhere. 



In such manner, the mean-field part (real-space pairing) 
and the fluctuation part ( pairing in k space) can be in- 
corporated thus into a single scheme. 
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